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ABSTRACT 

We investigate effects of aspherical energy deposition in core-collapse supernovae on the light curve 
of the supernova shock breakout. We performed two-dimensional hydrodynamical calculations of an 
aspherical supernova explosion to obtain the time when a shock wave generated in the stellar interior 
reaches the stellar surface in each radial direction. Using results of the calculations, light curves during 
the shock breakout are derived in an approximate way. We show that the light curve during the shock 
breakout can be a strong indicator of aspherical properties of core-collapse supernovae. 

Subject headings: X-rays: bursts — shock waves — supernovae: general 
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1. INTRODUCTION 

Massive stars with their main-sequence mass greater 
than 8M0 end their lives as core-collapse supernovae 
(CCSNe). The very beginning of dynamical evolution of 
a CCSN, at which a strong shock wave resulting from the 
collapse of the iron core abruptly emerges from the stellar 
atmosphere, is called the supern ova shock breakout. Ear- 
lier theore tical investigations ([Klein fc Chevalieil 119781 : 
[FaigiTOTl predicted that UV/soft X -ray photons hav- 
ing been produced in the stellar interior begin to es- 
cape from the stellar surface at the moment of the shock 
emergence. The enormous amount of the escaping pho- 
tons makes the moment the brightest phase throughout 
the evolution of a CCSN. Although observations of the 
shock breakout are very challenging due to its short du- 
ration (ranging from a few to a thousand seconds) , some 
cases of the detection of the shock b reakout have been 
reported, SN 2008D (ISoderberg et al. 2008; lUI [200l 

tLzzah e t all [2001 iModiaz et al.ll2q09'). S NLS-04D2dc 
:hawins kreral.l 120081: IGezari et al.l l2008h. and SNLS- 
06Dljd (,Gezari et al.ll2008[ ). The growing examples offer 
us good opportunities for comparing the theory of the 
shock breakout with observations. 

The nature of the shock breakout has been in- 
vestigated intensively in the framework of semi- 
analytical approaches based on se lf-similar solutions 
Jlmshennik & Nadczhin 198 ^ iMatzner fc McKed 
[l999t [Suzuki & Shigcvama 2010aJ) and nu merical 
simulations ([Shigeyama et ah 19881: iWooslevI 
Blinnikov k Nadyozhinl 119911: IE nsman fc Burrowa 



Tominaga et al.l 120091 ) . For the 
of the coupling of radiation and 
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proper treatment 
m atter, radiation- 
iShigevama et all 



Iwdrodynamics codes are often used 
([1988) performed radiation-hydrodynamical calculations 
based on the flux-limited diffusi on appro ximation. 
[Blinnikov fc Nad vozhinI ([1991 1 and Ensman fc Burrows 
([1992f ) developed multi-group radiation-hydrodynamics 
codes for the modeling of the phenomenon. Although 
there are differences between methods to model the 
phenomenon, these works reached the following general 
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agreements: i) the temperature behind the shock wave 
reaches about 10^ - 10® K, ii) the matter behind the 
shock wave emits UV/soft X-ray photons having a black 
body spectrum, iii) the light travel time across the 
progenitor star is crucial to determine the light curves 
of the emissions. It is noteworthy that the observed 
spectrum of the shock breakout of S N 2008D is well 
fitted by a power-l a w (ISoderber g ct al. 2008: Li '20081; 
iMazzah et all 120081: IModiaz et al...2009.) , which contra- 
dicts the theoretical expectation. Recent inves tigations 
([Wang et all 120071: ISuzuki fc Shigevamal [2010al) suggest 
that interactions between thermal photons and electrons 
at the shock front via Compton scattering , which is 
so-call ed bulk comptonization (Bla ndford fc Pavnd 
Il981allbh . can form a power-law spectrum. 

These investigations assume spherical symmetry, i.e., 
the shock breakout is assumed to occur simultaneously 
at every point on the surface of the progenitor star. How- 
ever, recent studies on the explosion mechanism of CC- 
SNe strongl y indicate aspherical deposition of the explo- 
sion energy (iMacFadven fc Wooslevlll999HBlondin et al.l 
[200aiKotake et al.ll2004ir although details of the energy 
deposition process remain unclear. From the observa- 
tional viewpoint, CCSNe are consid ered to be aspherical 
in general fe.g.. iMaeda et all 120081 ). For an aspherical 
explosion, the shock breakout does not occur simultane- 
ously in different radial directions. As a result, the light 
curve during the shock breakout of an aspherical explo- 
sion can deviate from those of spherical one. In other 
words, the light curve during the shock breakout can be 
used as a pr o be fo r the explosion geometry of a CCSN. 
iCouch et all ([2009f ) performed a series of hydrodynami- 
cal simulations of jet-induced explosions of a red super- 
giant progenitor. They investigated only the difference 
of light curves during the shock breakout between ther- 
mal and kinetic energy-dominated jets. However, light 
curves during the shock breakout must contain more in- 
formation on the explosion geometry of a CCSN, e.g., the 
viewing angle, the degree of aspherical energy deposition, 
and so on. 

In this letter, we present an approximate method to 
calculate light curves during the shock breakout assum- 
ing axisymmetric energy deposition and show that the 
thus calculated light curves indeed reflect the explosion 
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Fig. 1. — The emergence time as a function of the angle 9 for 
various a. The circles, squares and triangles are correspond to 
models with a = 0.8, 0.5, and 0.2, respectively. The solid, dashed, 
and dotted curves are calculated by the empirical formula J2j for 
a = 0.8,0.5, and 0.2. 



geometry of a CCSN. In §2, we describe the procedure to 
calculate light curves during the shock breakout. Resul- 
tant light curves are shown in §3. Finally, §4 concludes 
this paper. 

2. FORMULATION 

In this section, we describe our hydrodynamical models 
and the procedure to calculate light curves during the 
shock breakout. 

2.1. Hydrodynamical model 

As we noted in the previous section, the time when 
the shock wave propagating in different radial directions 
reaches the stellar surface (referred to as the emergence 
time in this paper) can take various values. To calcu- 
late light curves during the shock breakout, we must 
know the emergence time at each point on the stellar 
surface. Then, we performed a series of hydrodynami- 
cal calculations of a supernova explosion of a blue su- 
pergiant progenitor using a two-dimensional relativistic 
Eulerian h ydrodynamics code, which is d eveloped ac- 
cording to iDel Zanna fc BucciantTnil (|2002[ ). Details of 
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Fig. 2. — A schematic view of the geometry of the aspherical 
shock breakout. 



the calculations are described i n the forthcoming pa- 
per (jSuzuki fc Shigevamall2010"bl) . In the following, we 
present the brief summary of the calculations. The code 
integrates the basic equations of the special relativistic 
hydrodynamics in two-dimensional spherical polar coor- 
dinates (r, 9) , where r is the distance from the origin and 
9 is the angle measured from the symmetry axis. The 
equation of state for an ideal gas of the adiabatic index 
4/3 is a.ssumed. The progeni t or mo del is the same as 
used in iShigeyama fc Nomotol (|1990D . i.e., a progenitor 
model of SN 1987A. To generate a strong shock wave 
propagating in the stellar interior, we impose the follow- 
ing initial condition on the radial component of the 
velocity at the inner boundary (r = rin ~ 10^ cui), 



2-Bc: 



1 + q:cos(26I) 



47rr4Ari„pin ^1 - 2a/3 + 7a2/15 ' 



(1) 



where i^oxp = 10^^ erg is the explosion energy, is the 
density at the inner boundary, and Arin is the width of 
the radial zone at the inner boundary. Here we have 
introduced a parameter a characterizing the extent of 
the aspherical explosion. Using this parameter, we can 
express the ratio of the radial component of the veloc- 
ity at (rin,0) to that at (rin, 7 r/2) a s ?;in(0)/i;in(7r/2) = 
(1 + a)/(l - a). iMaeda et all (|2002D calculated the den- 
sity and the velocity distribution of ejected matter as 
a result of aspherical supernova explosions by a series 
of two-dimensional hydrodynamical calculations and ob- 
tained the distribution of h eavy element s , such as ^^Ni, 
synthesized in the ejecta. IMaeda et al.l ()2006[) studied 
the optical emission from the ejecta by using a Montc- 
Carlo radiative transfer calculation. They found that 
the axis ratio of the Ni blob in the ejecta is at most a 
few even for the model with Win(0)/win(7r/2) = 8, which 
corresponds to a = 7/9 in our model. Therefore, all 
models considered here must be observed as explosions 
with the axis ratio smaller than a few. It should be noted 
that there are many other candidates for the form of the 
angular dependence of the initial condition for t;i„. In 
this study, we adopt one of the simplest form, eq. (HJ, 
because it has an advantage in clarifying general behav- 
iors of light curves during the shock breakout originated 
from an aspherical explosion. However, due to this sim- 
ple treatment of asphericity, our model fails to repro- 
duce gamma-ray burst associated SNe even if we choose 
large a. Results of s imulations of jet-driven explosions 
fe.g.. lTominagal[2009[ ) show that matter located near the 
equatorial plane fall back to the central object, which 
strongly affects the velocity of the shock wave propagat- 
ing along the direction. To calculate the propagation 
velocity of the shock accurately, we need to perform sim- 
ulations of the explosion including self-gravity. For the 
boundary condition at the stellar surface, we assume the 
free boundary condition, i.e., the derivatives of physical 
quantities with respect to r at the stellar surface are as- 
sumed to be zero. 

The thus generated shock wave propagates outward 
and finally reaches the surface. Figure [T] shows the emer- 
gence time T{a, 9) as a function of the angle 9 for various 
a. The emergence time is defined so that T(a,0) — 0. 
The circles, squares and triangles in Figure [1] correspond 
to models with a = 0.8,0.5, and 0.2, respectively. We 
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find that the behavior of the emergence time in the range 
of < a < 0.8 is well approximated by the following em- 
pirical formula, 



T(a, 



65a(l + 3.6a) (sin^ 6 - 0.2 sin'^ 



(2) 



The lines in Figure [T] are plotted using this formula. 

To model the shock breakout accurately, we must treat 
the coupling of radiation and matter by solving the equa- 
tions of radiation-hydrodynamics, because the optical 
depth of the shock falls below unity after the shock emer- 
gence. However, thermal photons having been emitted 
by the shock front diffuse out from the envelope when 
the optical depth Ts of the shock decreases to the speed 
of light divided by the shock velocity, c/Vs- In the case 
of the explosion of a blue supergiant progenitor with the 
explosion energy of i?exp = 10^^ erg, the shock velocity 
reaches ~ 0.03c at the moment of the shock emergence, 
which leads to the optical depth at the shock breakout 
of Tg ~ 30. The radiation field at such a large optical 
depth is expected to approach the black body. Further- 
more, in the Rosseland approximation, the fraction of 
the anisotropic part of the radiation field is of the order 
of 4/(37rTs), which is much smaller than unity. Thus, 
we assume that isotropic black body radiation with a 
constant flux lasts for a time interval r after the shock 
reaches the stellar surface. We need the temperature 
of the radiation emitted by the matter behind the shock 
wave when the shock reaches the surface in order to eval- 
uate the flux at each point on the stellar surface. In 
the hydrodynamical calculations, we assumed the equa- 
tion of state for an ideal gas of the adiabatic index 4/3, 
which is a good approximation if the pressure is radia- 
tion dominated. Therefore, we estimate Tj- by assuming 
that the post-shock pressure Ps at the outer boundary 
is equal to the radiation pressure given by where 
a(= 7.566 X 10^^^ erg cm^"^ K^**) is the radiation con- 
stant. We flnd that the post-shock pressure when the 
shock reaches the stellar surface is almost independent 
of the angle 9, Ps — 2.4 x 10^ g/cm s^, which leads to the 
radiation temperature of ~ 10^ K. 

2.2. Derivation of light curves 

We calculate light curves during the shock breakout for 
a given emergence time distribution T(a, 9). A schematic 
view of the geometry considered here is shown in Figure 

m 

At flrst, we introduce the viewing angle O as the 
angle between the line of sight and the symmetry 
axis. We define cartesian coordinates {x,y,z) so that 
the z-axis is identical with the line of sight. The 
symmetry axis is assumed to be in the x-z plane. 
Thus, a unit vector parallel to the symmetry axis 
is expressed as (sin 0, 0, cos 8) in the coordinate sys- 
tem. We consider the flux F{t,9o,(l)o,a) at a point 
{R sin ^0 cos 0o i R sin sin c^o j R cos ) on the surface of 
a progenitor star with the radius 3 x 10^^ cm) at 
a time t. Using the angle x between the symmetry axis 
and the position vector of the point given by 



cos X = cos ^0 cos 6 + sin sin Q cos (j)o , 



(3) 



we can express the emergence time at the point as 
r(a, x)- As we noted in the previous section, we assume 
that black body radiation with a constant flux lasts for 



a time interval r after the shock emergence. The time 
interval r is the diffusion time scale of photons emitted 
from the shock front. The value depends on the explosion 
energy, the mass of the ejected matte r ^ and the structure 
of the envelope. 'Matz ner fc McKed (1999) studied the 
shock brea kout emission b y using the self-similar solution 
derived bv lSakurail (|1960f ) and conflrmed that the results 
are consistent with those of numerical simulations. They 
derived 40 s for the explosion of a blue supergiant pro- 
genitor with the explosion energy of 10^^ ergs and the 
ejecta mass of IOM0. Thus, we set r = 40 s in this 
study. The flux is expressed as 



F{t, 90,(1)0, a) 







for T{a,x)<t<Tia,x) 
otherwise, 

(4) 

^ is the Stcfan- 



T, 



^ erg cm ^ s 



K- 



where ctsb = 5.67 x 10^ 
Boltzmann constant. 

Next, we consider an observer located at (0,0,1?). 
The distance D is sufficiently larger than the ra- 
dius R, i.e., the progenitor star is regarded as a 
point source for the observer. A photon observed 
at a time t was emitted from the source position 
(i?sin^oCOS0o,-Rsin0osin(/)o,i?cos6'o) a.t to = t — {D — 
Rcos9o) /c. Therefore, the luminosity L{t,a,Q) at the 
position of the observer at t is obtained by integrating 
the flux F{to, 9o,(t>o) multiplied by R^ sin^o with respect 
to < 6*0 < 7r/2 and < (j}o < 27r, 



L{t,a,Q) = R^ 



7r/2 



F{tQ,9o,(f>Q,a) sin 9od9od(t)o- 



When we approximate the form of the emergence time 



as T{a, 



A{a) sin 9 and assume = 0, we can carry 



out the integration in the above formula ([5|) and derive 
the following approximate formulae of the luminosity. 



L(i,a,0) 



2^R^asBT^[l- f{a,t)] 
for ^ <t<i 

2TTR^asBTf[fia,t-T) 
for - 

2TTR^asBT,^fia,t 



T<t<^ 

-t) 







for f + A{a) <t<^ 
otherwise. 



where we have defined a function f{a,t) as 



r, 
t)] 

f yl(a), 
T + A{a), 
(6) 



/(a,i) 



-R + ^/R^ - 4cA(a)[ct - cA{a) - D] 
2cA{a) ' 



(7) 



We find that the above approximate formula agrees with 
numerically calculated light curves within 10 percent 
error if we define the function A{oi) as 



A(a) = 52a(l -1-3.60;). 



3. RESULTS 



(8) 



We calculate light curves during the shock breakout by 
using cqs. ©, ©, (HI), and dSJ. The distance D is set 
t o be £> = IQQR. 

lEnsman fc Burrows! ()1992D . who tackled the problem 
with a multi-group radiation-hydrodynamics code, cal- 
culated light curves of the breakout emission for a similar 
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Fig. 3. — Light curves during the shock breakout for various 
viewing angles © and a fixed a(= 0.5). The solid, dashed, dotted 
and dash-dotted lines represent models with = 0°, 30° , 60° , and 
90°, respectively. 
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Fig. 4. — Light curves during the shock breakout for a fixed 
viewing angle © = 0° and various ct. The solid, dashed, dotted, 
and dash-dotted lines represent models with ol = 0.8,0.5,0.2, and 
0, respectively. 
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Fig. 5.— Same as Figure H] but for 6 = 90°. 



reaches the peak value and then begins to dechne. 

This behavior is easily explained by considering how 
photons emitted from the stellar surface travel to the 
observer. In the case of small viewing angles, the shock 
wave initially emerges from the point on the symmetry 
axis {9o ~ 0) and then points with > experience the 
shock emergence one after another. Because the emer- 
gence time distributions T{a,6) flattens at 9 = 7r/2 as 
shown in Figure [H the shock emergence at points with 
9o ~ 7r/2 is expected to occur almost simultaneously. In 
other words, in a later epoch, the observer receives more 
photons than in the earlier epoch. This is the reason why 
the luminosity reaches the peak value after a gradual in- 
crease for small viewing angle. 

Next, we consider models with large viewing angles. In 
this case, the shock wave initially emerges from points 
with 6*0 — tt/2. Photons emitted from the points need 
longer time to reach the position of the observer than 
those from points with 9o < tt/2. Because of the com- 
bined effect of the time delay and the early emergence 
of the shock, photons from points with 9q c± tt/2 reach 
the position of the observer almost simultaneously, which 
leads to a more sudden brightening to the peak luminos- 
ity. 



progenitor model. Some general features of their results 
are: i) the peak luminosity reaches ^ 10^^ erg/s, ii) the 
duration of the emission is a few minutes, iii) the radi- 
ation temperature is of the order of 10^ K. Despite the 
approximate treatment of emissions from the shock, light 
curves calculated by our method successfully reproduce 
the features. 

3.1. The dependence on the viewing angle 

Figure[3]shows light curves for various viewing angles Q 
and a fixed a{= 0.5). A clear distinction between models 
with small viewing angles and those with large viewing 
angles is recognized. For small viewing angle models, 
the luminosity gradually rises after a sudden brightening 
and eventually reaches the peak value. For large viewing 
angle models, on the other hand, the luminosity suddenly 



3.2. The dependence on the parameter a 

Next, we investigate how the degree of asphericity of an 
explosion, which is characterized by the parameter a in 
eq. ([IJ, affects the light curve during the shock breakout. 

Figure |4] shows light curves for a fixed viewing angle 
8(= 0°) and various a. The dash-dotted line represents 
the model with a = 0, i.e., a spherical explosion. In this 
case, the shock breakout simultaneously occurs in the 
same way at every points on the stellar surface. There- 
fore, the luminosity does not evolve with time. The du- 
ration of the emission is roughly determined by light- 
travel-time across the stellar radius, R/c ~ 100 s. For 
models with a larger a, the duration of the emission is 
longer than that of the model with a — 0. This is be- 
cause the shock emerges at points with 9o > after it 
does at the point with 6*0 = 0. Therefore, the duration 
of the emission is roughly estimated as R/c + T{a,TT/2). 
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Furthermore, since the total energy of the burst is same 
for all models, the average luminosity is small for models 
with large a. 

Figure [5] shows light curves for a fixed viewing angle 
8(= 90°) and various a. In this case, the duration of 
the emission takes a similar value for all models. For 
models with a large a, the feature claimed above, i.e., 
the sudden rise and the subsequent gradual decrease of 
the luminosity, is prominent. 

4. CONCLUSIONS 

In this letter, we investigate effects of an aspherical ex- 
plosion on the light curves during the supernova shock 
breakout. To obtain the emergence time distribution 
T(a,9), we performed a series of two-dimensional hy- 
drodynamical calculations of aspherical supernova explo- 
sions. We establish an approximate way to calculate light 
curves during the shock breakout using results of the sim- 
ulations and confirm that resultant light curves success- 
fully reproduce some features of the emission from the 
supernova shock breakout. 

We find that the light curves strongly depend on the 
viewing angle and the degree of asphericity of an explo- 
sion. For small viewing angles from the symmetry axis, 
the luminosity gradually increases and then reaches the 
peak value. On the contrary, for large viewing angles, 
the luminosity suddenly reaches the peak value and then 
gradually decreases. For a larger degree of asphericity, 
the light curve can be distinguished from that of a spheri- 
cal symmetric explosion irrespective of the viewing angle. 
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These features allow us to use light curves during the su- 
pernova shock breakout as a probe for the explosion ge- 
ometry of CCSNe. It is noteworth y that the light cu rve of 
the shock breakout of SN 2008D (ISoderberg etaI1[200l 
is similar to that of models with large viewing angles, 
e.g., the model with a — 0.5 and Q — 90, except for the 
duration. The nebu lar phase observations of SN 2008D 
(|Tanaka et al.ll2009[ ) revealed that SN 2008D was a side- 
viewed bipolar explosion, which also prefers large viewing 
angle models. 

In this study, we only calculate the supernova shock 
breakout from the progenitor of SN 1987A, i.e., a blue 
supergiant, in order to confirm the validity of our ap- 
proximate way to calculate light curves by comparing re- 
sults of other studies that adopt more sophisticated ways. 
However, since the duration of the shock breakout from 
red supergiant progenitors, whose radii are typically of 
the order of 10^"^ cm, is expected to b e longer than those 
from blue supergiant progenitors (jMatzner fc McKei 
|1999( ). the shock breakout from a red supergiant pro- 
genitor is likely to be obs erved by future observations . 
In the forthcoming paper (|Suzuki fc Shigevamal l2010b'). 
we will present results of the calculations for other pro- 
genitor models as well as the description of details of the 
hydrodynamical simulations. 
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